Prom avalanches to fluid flow: A continuous picture of grain dynamics down a heap 
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Surface flows are excited by steadily adding spherical glass beads to the top of a heap. To simul- 
taneously characterize the fast single-grain dynamics and the much slower collective intermittency 
of the flow, we extend photon-correlation spectroscopy via fourth-order temporal correlations in the 
scattered light intensity. We find that microscopic grain dynamics during an avalanche are similar 
to those in the continuous flow just above the transition. We also find that there is a minimum 
jamming time, even arbitrarily close to the transition. 
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We have all poured sand in a heap. If grains are added 
rapidly, the heap relaxes through continuous, fluid- like 
flow at its surface. If grains are added slowly, it re- 
laxes intermittently through discrete avalanches. The 
existence of both continuous and intermittent regimes 
is characteristic of the dual nature of granular media, 
which can exhibit properties of both liquids and solids 
Jl],Q. Similar behavior, such as bubble rearrangements 
in foams or kinetic heterogeneities in supercooled liquids, 
arises in other systems near the onset of jamming ||Q. 
In spite of its ubiquity, and its obvious importance for 
applications and geophysical phenomena, the transition 
between intermittent and continuous flow has received 
relatively little attention. One reason may be that ki- 
netic theories ||,|J of granular hydrodynamics, as well as 
minimal- ingredient theories of surface behavior [JtJ , apply 
only to continuous flows, while cellular-automata models 
of avalanches apply only to rate- independent intermittent 
flows j^]; both approaches break down near the transi- 
tion. Another reason may be the vast separation in time 
scales between the grain velocity fluctuations and the in- 
termittency dynamics, which makes it difficult to capture 
both phenomena simultaneously. 

Recently we used diffusing- wave spectroscopy (DWS), 
a photon-correlation spectroscopy for opaque media J9|, 
to probe motion in homogeneously excited sand p0|. 
DWS reveals dynamics on time- and length-scales too 
small for imaging, but cannot be applied to intermittent 
flows. This is because in the static periods all the grains 
are perfectly correlated, violating the assumption that 
the total scattered electric field has Gaussian statistics. 
Here, for a heap upon which grains are steadily poured, 
we show how intermittency can be detected and char- 
acterized using higher-order temporal correlations in the 
intensity of multiply-scattered laser light. We extract 
single-grain dynamics, as well as probability functions 
describing the collective intermittency of the flow. This 
gives an unprecedented picture of the dynamics in both 
intermittent and continuous flows, as well as across the 
transition between these extremes. 

System. The experimental setup (Fig. Qa) consists of 
parallel static-dissipating plastic walls, 30x30 cm 2 , form- 



ing a w = 9.5 mm wide vertical channel closed at the 
bottom and on one side. We use spherical glass beads of 
diameter 0.33±0.03 mm, large enough to minimize effects 
of interstitial gas and electrostatic forces while keeping 
the system size manageable; the density of the packing 
is p = 1.35 g/cc. Flow is created by dropping grains into 
the closed side of the channel, and breaking their fall by 
a secondary hopper just above the heap. Grains falling 
out the open end are recycled by a vacuum-driven ele- 
vator. The flow rate is set by a micrometer-controlled 
gate valve and a knife-edge, allowing for wide variation 
of 0.025-5 g/s (50-10,000 grains/s). Prior to a run, the 
beads are washed, dried in air, and circulated until a 
steady state is reached. 
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FIG. 1. The experimental setup (a), and two macroscopic 
measures of the response vs mass flow rate: (b) heap angle 
and (c) layer thickness. 
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For orientation the macroscopic behavior of the heap 
is summarized in Figs. [l]b-c vs mass flow rate, Q, the 
sole parameter that sets the grain dynamics. As opposed 
to the inclined-plane geometry (IT), our system chooses 
both the angle 9 of the heap and the thickness 5 of the 
flowing layer. This avoids such issues as the role of the 
plane roughness, and should ultimately permit more di- 
rect contact with theory. Visual inspection reveals that 
the transition between intermittent and continuous flows 
is at Q c = 0.36 g/s, and that it is rather sharp. By con- 
trast with the rotating-drum geometry (f^,[l3| , this tran- 
sition is not hysteretic |l3]] , and avalanches always start 
at the top. Also by contrast, both 5 and 9 are uniform 
across the whole surface. Fig. |f| shows that these quanti- 
ties become independent of Q in the intermittent regime, 
with 9 ranging between a minimum angle of repose and 
a maximum angle of stability, and with S being the same 
thickness as at the transition, 5 C — 5 mm, which is about 
15 grain diameters or about half the width of the chan- 
nel. Thinner layers cannot flow indefinitely, as dissipative 
forces overwhelm gravity and inertia. Fig. [l] also shows 
how both 9 and S increase monotonically for continuous 
flows above the transition. Interestingly, the heap angle 
varies logarithmically with Q and, as measured at the 
walls for Q = 0.5 g/s, the velocity vs depth is nearly ex- 
ponential, v x (z) « (1 cm/s) exp[— z/(0.15 cm)] for z < 5 
This is unlike the modified-Bagnold mod el p redic- 
tions of 6» - 6» c oc Q 2 and v x (z) oc [5 3 ' 2 - z 3 / 2 ) pf. This 
implies that there must be enduring contacts and/or that 
the rate of collisions is not set by dv x /dz. 

Light Scattering. While visual measurements provide 
a rough overview of a heap's macroscopic response, 
multiple-light scattering measurements provide a de- 
tailed picture of both individual and collective grain dy- 
namics over a broad range of time scales. Here we illumi- 
nate the top surface with an Ar + laser, wavelength A = 
514 nm, over a 8 mm region located L = 18 cm from the 
closed end. The intensity fluctuations in the speckle pat- 
tern created by backscattered photons are captured by a 
single-mode optical fiber coupled to a photon-counting 
assembly. Following Ref. |15j , we measure simultane- 
ously in real-time both the second- and fourth-order in- 
tensity correlation functions g^ 2 \r) = (J(O)J(r)) / (I) 2 

and g^(r) = (I(0)I(T)I(t)I(t + T)> / (if, where T is 
a fixed delay that can be chosen between 50 ns and 50 ms, 
and r is swept over a range of values by a digital correla- 
tor. In analogy with the 'second-spectrum' |Q, g^\r) 
can be thought of as a 'second-correlation' or 'correlation 
of correlation' quantifying fluctuations in the correlation 
signal I(t)I(t + T) with time t. 

Examples of g^ 2 \r) and g^\r) are shown in Figj^. 
For continuous flows, both correlations exhibit a single 
complete decay at short time scales. Good agreement is 
seen between (t) data and the predictions generated 
solely from g^ir) data, as in Ref. |Lj|, showing that 
the fluctuations are Gaussian. Therefore, the usual ma- 
chinery of DWS may be applied to extract single-grain 



dynamics. For intermittent flows, by contrast, both cor- 
relations exhibit a partial decay at short times, followed 
by a plateau and a final decay at long times. Further- 
more, the fluctuations are non-Gaussian, implying that 
a new model is required to extract grain motion. 
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FIG. 2. Examples of (a) second- and (b) fourth-order 
temporal correlations in the intensity of laser light multi- 
ply-scattered from the surface of the heap. For continuous 
flows, the fluctuations are Gaussian and the correlations ex- 
hibit a single decay; for intermittent flows, the fluctuations are 
non-Gaussian and the correlations exhibit a two-step decay. 

Several ingredients are required to model the full spec- 
trum of grain dynamics. One is the usual normalized 
electric field autocorrelation function 7(7"), whose Fourier 
transform gives the broadening of the power spectrum 
due to relative motion of scattering sites during flow. 
For our backscattering geometry, the theory of DWS 
|,0 gives | 7 (r)| oc exp{-2^k 2 6v 2 T 2 + W/Q, where 
k = 2tt/X, Sv is a typical relative grain speed averaged 
over the scattering volume, I* = 1.3 mm is the photon 
transport mean- free path, and l a = 1.0 cm is the pho- 
ton absorption length. Since Sv is the only quantity not 
known in advance, it may be extracted from the data. 

All other ingredients relate to the probabilities with 
which the system switches between flowing / '1' and 
static / '0' states. Most important is the probability 
Po( r ) f° r the system to be in the same static config- 
uration after a time interval r. This function decays 
monotonically from one to zero. It can be related to the 
probability P(r) that two flowing states are separated 
by r via Pq(t) — P(t')(t' — T)dr'/t , an d hence by 
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P(t) — t d 2 Po/dT 2 , where t is the average static state 
duration. Also important are the probabilities Pij(r) for 
the system, if initially in state i, to be in state j after a 
time t. These interrelated functions decay from either 1 
or to the fraction of time fj spent in state j. 

Making the simplest assumption consistent with visual 
observations, that starting/stopping transients arc negli- 
bible and that the speckle pattern fluctuates rapidly dur- 
ing each flowing state but is constant during each static 
state, we find 

fl C0(r) = 1 + /3 [hP 11 {r)\ 1 {r)\ 2 + f P (r)} , (1) 

where 1//3 is roughly the number of speckles viewed by 
the detector. For continuous flow, where f\ and Pn are 
identically one, Eq.(j|) reduces to the usual Siegert rela- 
tion connecting intensity and field autocorrelations. For 
intermittent flow, Eq.(|lj) exhibits a decay at short times 
due to j(t), during which the switching functions are 
constant, and a final decay at long times due to Pq(t). 
This is unlike the two-step decay seen in supercooled sys- 
tems, where the field fluctuations are Gaussian and the 
final decay (the a-relaxation) is due to cage-hopping ef- 
fects. Note that the intercept g^ 2 \0) = 1-1-/3 gives the 
same second intensity moment as for a Gaussian process; 
this property also holds in our general expressions for 
and g( 4 \ Intermittency gives non-Gaussian fluctuations 
because of grain correlations during each static state, 
unlike the broader intensity distributions familiar from 
number fluctuations in single-scattering experiments. 

The predictions for higher-order correlations are cum- 
bersome, but reduce considerably for the usual case of 
good separation in decay scales for j(t) and the switch- 
ing functions. Choosing an intermediate value for the 
fixed delay T, so that 7 (T) = and P(t ± T) = P(r), 
we find 

^ 4) (t)«1+/3/i[2 + /?|7(t)| 2 ] \l{r)\ 2 + 

/3 f [2 + pP 00 (t) + (4 + 10/3 + 6/3 2 ) P (t)] . (2) 

This neglects terms in j(\t — T\), which cause a peak 
in g}^ \t) at t = T that is characteristic of our slice 
of the full g^(ri, t 2 , t 3 ). As above, Eq.(||) exhibits a 
single decay for continuous flows and a two-step decay 
for intermittent flows. Armed with Eqs. (|l|j^), one can 
not only reliably detect the presence of intermittency, 
but can also characterize both the grain motion via 7(7-) 
and the collective intermittency of the flows via Po(t) 
and Pjj(r). This extends the utility of photon-correlation 
techniques (not just DWS) to systems like slowly-driven 
sand exhibiting intermittency. 

Grain Dynamics. For continuous flows, according to 
the above discussion, the correlations fully decay at a 
time scale given approximately by X/Sv. Using Sv as 
the sole parameter, we obtain excellent fits to g^(r) 

(and equivalently to grp (r) since here the fluctuations are 
Gaussian) as shown in Fig. ||(a). The results in Fig. |^(a) 
show that Sv varies logarithmically with the flow rate, 



just like the heap angle. The value of Sv appears to re- 
flect both shear flow and random ballistic motion. At 
Q = 0.5 g/s, slightly above the transition, the r.m.s. av- 
erage strain rate in the scattering volume is estimated 
from v x (z) and Q to be less than 4/s. This gives a con- 
tribution to Sv of l*y/(j 2 ) /5 = 0.4 cm/s Q, which ac- 
counts for half the fitted value of 0.8 cm/s. We believe 
that the random contribution is not due to free-flight 
between grain-grain collisions, since the Bagnold model 
does not apply, but is rather due to jostling parallel and 
perpendicular to the average flow as shearing layers slide 
over one another with grains always maintaining a near 
close-packing arrangement. While there is no saltation, 
altered dynamics at the free surface may also contribute. 
In any case, the extracted values of Sv represent a mea- 
sure of microscopic velocity fluctuations for the flowing 
grains. These are related to the average flow speed as 
Sv ~ w o.5±o.i^ smi ii a r to previous observations for a ver- 
tical hopper 0. The increase of fluctuations relative to 
the average may be the ultimate cause of jamming at low 
driving rates. 
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FIG. 3. Grain dynamics vs mass flow rate: (a) velocity 
fluctuations, (b) fraction of time spent static, (c) average time 
spend static, and (d) average time spend flowing. Consistent 
results are obtained from second- and fourth-order correlation 
data using Eqsjl] and [| respectively. 

For intermittent flows, according to the light scatter- 
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ing discussion, the correlations partially decay at X/6v 
and then fully decay at a time scale reflecting the inter- 
mittency dynamics. Eqs. (^|]|) show that between these 
decays the correlations plateau at values of 1 + /3fo and 
1 + /3(6 + 11/3 4- 6/3 2 )/o, respectively. Knowing f3 « 1, we 
easily extract the fraction of time fo the system spends 
at rest. The results in Fig. ||(b) show that the g^ 2 \r) and 
glp (t) plateaus are consistent and in accord with expec- 
tations. As the flow rate goes to (Q c ), fo approaches 1 
(0), denoting a completely static (fluid) state. 

Given the correlation plateaus, we may also fit the 
early-time decay for grain-scale velocity fluctuations dur- 
ing avalanches using Eqs.([j]-^). The results in Fig. ||(a) 
show that 8v is approximately independent of Q in the 
intermittent regime. Furthermore, the value is close to, 
but perhaps slightly less than, that for fluctuations in 
continuous flow just above the transition. So avalanches 
are not very different from barely-continuous flow. 

Finally we analyze the late-time decay of the inten- 
sity correlations in terms of avalanche statistics. Specif- 
ically, the functional form of the the switching function 
Po(t) may be deduced from g^(r) data using Eq.(jl]) 
and the known plateau values. Similarly Pq(t) may be 
deduced from (r) data, ignoring the small contri- 
bution in Eq.(B) from Pqq(t). The two results agree, 
and for all Q are consistent with the empirical form 
P ( T ) = exp(-r/t )/[l + (T/t ) 2 /2\. Namely, the first 
half of the decay is linear, and the second half is faster 
than exponential. This rules out self-organized critical- 
ity, for which Po(t) would be a power-law, as well as 
random-telegraph switching, for which Pq(t) would be 
perfectly exponential. Rather, the distribution of inter- 
vals between avalanches is peaked around an average, to, 
similar to behavior seen in Ref. [fiJJ. 

Results for the average static and flowing times, de- 
duced respectively from Po(t) = 1 — r/io + ... and 
fo = ta/(to + ti), ar e shown in Fig. ||(c-d) as a func- 
tion of flow rate. For the smallest Q, we find evidence of 
a quasi-static limit in which the time between avalanches 
scales as 1/Q and in which the duration of avalanches 
t\ approaches a constant. The observed times are on 
the order of the upper limits pwLS c /Q and pS c wL/Q c 
needed to fill and to empty, respectively, a volume wLS c 
of sand above the illumination spot. As the flow rate is 
increased toward Q c , where continuous flow commences, 
the time spent flowing appears from Fig. |^(d) to diverge. 
Interestingly, at the same point, the time spent in the 
static state does not vanish continuously but rather ap- 
proaches a nonzero constant. Even arbitrarily close to 
smooth flow, it takes significant time for the system to 
start moving again after it has jammed. This is analo- 
gous to common experience on a freeway, where traffic 
jams develop very suddenly but clear much more slowly. 

Conclusions. In summary we have introduced a novel 
optical method and used it to probe, simultaneously, 
both avalanche statistics and microscopic grain dynamics 
within the avalanches. A unique feature is that the tech- 



nique applies to both continuous and quasi-static inter- 
mittent flows, as well as to the transition between these 
extremes. This should enable a wide range of new ex- 
periments on granular dynamics as a high driving rate of 
tumbling, shaking, blowing, sprinkling, shearing, etc., is 
decreased and the system develops intermittent dynamics 
as a precursor to complete jamming. For example we are 
currently investigating the influence of channel dimen- 
sions, and the behavior as a function of depth, for the 
sprinkling geometry presented here. The results should 
help guide and test fundamental theories of granular flow 
that are based on the actual micrsocopic grain dynamics. 
This work was supported by NSF grant DMR-0070329. 
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